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SUMMARY 

An examination of the effects of compressibility, variable proper- 
ties, and body forces on fully developed laminar flows has indicated 
several limitations on such streams. 

In the absence of a pressure gradient, but presence of a body force 
(e.g. , gravity), an exact fully developed gas flow results. For a liquid 
this follows also for the case of a constant streamwise pressure gradient. 
These motions are exact in the sense of a Couette flow. In the liquid 
case two solutions (not a new result) can occur for the same boundary 
conditions. An approximate analytic solution was found which agrees 
closely with machine calculations. 

In the case of approximately exact flows, it turns out that for 
large temperature variations across the channel the effects of convec- 
tion (due to, say, a wall temperature gradient) and frictional heating 
must be negligible. In such a case the energy and momentum equations 
are separated, and the solutions are readily obtained. If the temper- 
ature variations are small, then both convection effects and frictional 
heating can consistently be considered. This case becomes the constant- 
property incompressible case (or quasi-incompressible case for free- 
con vection flows) considered by many authors. 

Finally, there is a brief discussion of cases wherein streamwise 
variations of all quantities are allowed but only in such form that the 
independent variables are separable. For the case where the streamwise 
velocity varies inversely as the square root of distance along the chan- 
nel, a solution is given. 
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INTRODUCTION 

Among all possible fluid flows, one of the most useful is the fully 
developed (i.e., independent of streamwlse distance) ch ann el flow. The 
flow is taken to he the motion generated by a constant pressure gradient 
(the familiar Poiseuille flow) or by a body force (ref. l) . In either 
case one us ually considers only an incompressible or quasi-incompressible 
flow with fixed fluid properties. This is in marked contrast to the case 
of Couette flow where two parallel walls move with respect to each other. 
In such a case there is no need to limit oneself to a perfect gas or to 
any particular variation of the transport properties (ref. 2 is a case 
in point) . 

The crucial difference between the Couette and the Poiseuille flows 
is that the former admits a stream wherein nothing depends on the stream- 
wise distance, while the latter requires that the pressure vary in the 
flow direction. Hence, to some small degree, at least, the other fluid 
properties will also vary in that direction if the state equation in- 
volves the pressure. The present study is an examination of the general 
circumstances under which there can be a fully developed laminar flow 
past fixed boundaries. 

One special problem considered is an unusual situation found by 
Ostrach (refs. 1 and 3 to 6) . He discusses the flow of a fluid in a two- 
dimensional channel under the influence of gravity. Incompressible flow 
is assumed except as is. required to generate a varying body force, anrl 
the fluid transport properties are assumed not to vary. Under these 
assumptions , the surprising result is found that there are two solutions 
to the flow in question for a certain range of values of the flow param- 
eters. The first corresponds roughly to the neglect of frictional heat- 
ing, while the other is near the (nontrivial) solution for homogeneous 
boundary conditions. 

There may be some doubt as to the stability of one of the solutions, 
presumably the second one. In any case, a question arises about the ef- 
fect of considering a real fluid having variable viscosity and ther mal 
conductivity as well as being truly compressible. 

Accordingly, the present paper treats the consequences of such gen- 
eralizations. However, to reiterate, one serious restriction is on 

all the flows considered herein: The flow is always fully developed, 

with the result that the effects of conditions near either end of the 
channel are ignored. 
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ANALYSIS 


Consider a two-dimensional flow of a viscous compressible fluid 
acting under the influence of an axial body force such as gravity. Var- 
iable viscosity and thermal conductivity are ad mi tted. The configura- 
tion is shown in the following sketch: 
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(pu) x + (pv) Y = 0 

p(ui^ + VUy) + P X = -Pf + 2G«fc) X + |^(uy + Vjr)^ - |Jp(^; + x 

P(UV X + Wy) + Py = 2(pVy)y + ^(Uy + V X^ X “ f + Y 

[ 222 
2Uy + Uy + V X + 

2 2 21 
2vy + 2uyV x - (ujr + v y) 


(1) 

( 2 ) 

(3) 


(4) 


(Symbols are defined in appendix A. ) Consider, in addition, two possible 
forms of a state equation, one applying for a gas, the other, for a 
liquid: 


P = pRT 

p = p [l - P(T - T)J 


(Gas) 

(Liquid) 


(5a) 

(5b) 
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In equation (5b)., p is the (small) volumetric expansion coeffi- 
cient , and p and T are reference values. The significant difference 
between the two state equations is that the second is independent of 
pressure. 

The boundary conditions on the channel walls are 

u(X,0) = u(X,d) = v(X,0) = v(X,d) = 0 \ 

T(X,0) = T^X) ( (6) 

T(X,d) = T^X) ] 

where d is the distance between the channel walls. The temperature 
boundary conditions could, of course, be replaced entirely or in part by 
a heat-transfer condition but, for the purpose of this report, such a 
change is unimportant. 

Equations (l) to (6) are sufficient to define the fully developed 
flow in a channel provided the viscosity and conductivity variations with 
temperature are known and provided further that the forced-flow pressure 
gradient, if any, is specified. In seeking solutions of these equations 
for flows in a very long channel (i.e., fully developed) we are led to 
consider three approaches: first, exact solutions entirely independent 

of distance (X) along the channel; second, solutions approximately inde- 
pendent of X; and third, solutions wherein the variables are separable. 
In each case the results can be expected to differ according to which of 
the state equations applies. In this connection, it is important to ob- 
serve that the viscosity and thermal conductivity vary differently in 
liquids and gases. In particular, the viscosity rises with temperature 
for a gas and falls for a liquid. 

The phrase "exact solution" should perhaps be defined. In this re- 
port it is understood to mean a solution which satisfies equations (l) 
to (6) rigorously. However, no consideration is given to conditions 
near the ends of the channel. There are two relatively simple circum- 
stances under which such exact solutions can be found. For a gas nothing 
can vary with X, not even the pressure. For a liquid this restriction 
is moderated to the extent that only the gradient of the pressure need 
be independent of X. This relaxed condition occurs because of the 
pressure-independent state equation (eq. (5b) ) for a liquid. 


Exact Flow of Gas 

If nothing depends on X, equations (l) to (5a) and (6) become 


v = 0 


(7a) 
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(m-u y )y = P f 
P y = 0 

(^T y ) y = -pu^ 
P = pRT 

u(0) = u(d) = 0 

*(°) = nv 0 

T(d) = 


(7b) 

(7c) 

(7d) 

(7e) 

( 8 ) 


Assume that the viscosity and thermal conductivity vary as powers 
of temperature, and also, for convenience, change the independent 
variable (Y). 


Thus , suppose 


and let 


or 



(9) 


( 10 ) 
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where the constant B is as yet undetermined. Then equations (7) 
become 





If the viscosity varies linearly with temperature (n-j_ = l) , then these 
equations are separated. If ng is also unity, the solution is 



Finally, B can be found from the second form of equation (lO) . Thus 


or 


where 
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(14) 


(14a) 
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The quantity in brackets is essentially the parameter K defined 
in reference 1. This is always positive. Under these circ ums tances 
equation (13) has only one real root, that root being such that 


0 < 


HD, 


m 


dT w, 


< 1 


For example, for air under standard conditions, if d = 3 and f is 
gravitational acceleration, equation (14) yields BT^/dT,^ = 0.904. A 

convenient standard for comparison of the present solution with more 
approximate results is given by the mass flow. This is 



(15) 


The corresponding result for incompressible flow with constant fluid 
properties, or for compressible flow with variable properties, but ne- 
glecting frictional heating, is 



For the cane cited prior to equation (15) , where ET^/dl^ = 0.904, the 

actual mass flow is about 25 percent below the incompressible value. 
Unfortunately, the mean flow velocity is about 800 feet per second. This 
high velocity would probably preclude the possibility of laminar flow 
even existing. If an example leading to slower flow were considered, the 
difference between the two results would have been small. This demon- 
strates that the variation of viscosity is unimp ortant. 


A final comment: If the viscosity is not assumed to vary linearly 

with the temperature, the momentum and energy equations cannot be sepa- 
rated. As this circumstance (lack of separation) led Ostrach (refs. 1 
n. nd 4 to 6) to find two solutions for the flow rather than one, it is 
perhaps worth examining further. Define 



( 16 ) 
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Then equations (ll) become 


u 




aB^Pf /aH(n2 ~ n l + 1 ) \ ^ n ' 1 ^^ 1 n l +n 2^ 

\ b / 


P = Constant 

H nn - "A 


(17a) 


(17b) 


For gases, one expects that n^_ < 1 and > 0; "thus the right side of 

equation (17a) is a decreasing function of H. In this case there can^ 
be at most one solution of the problems. The argument goes this way. 
Suppose one solution is known. If a second solution has larger H, then 
by equation (l7a), u^^ is smaller. Hence, for a reasonable profile, u^ 

is reduced. Then, by equation (17b), is of lesser magnitude. Hence 

for a reasonable case, H must also be small, which is a contradiction. 


It should be emphasized, however, that this case of no pressure gra- 
dient whatever has no connection with the work reported in references 1 
and 4 to 6. 


This completes the solution for the exact fully developed flow of a 
gas. There are two generalizations which can readily be made. These in- 
volve the addition of a body force transverse to the channel and the ad- 
dition of heat sources in the fluid. The solutions are given in appendix 
B. 


Exact Flow of Liquid 

Here it is assumed that there is a pressure gradient such that, at 
least, P X and Py are independent of X. No other X- dependence is ad- 
mitted. Then the system is that given in equations (7) and (8), except 
for the X-momentum equation, which can conveniently be written as 

(nuy)y = -pp(T - T)f + (F x + pf) - pPf(T - T*) (18) 


where 

* ?x + P(1 + 


P0f 


( 19 ) 
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It is -worth, observing here jbhat , if the reference point is changed in 
equation (5b)(i.e., a new T) , this has no effect whatever on. the 
value of T*. This is because the state equation (eq. (5b)) is really 
of the form p = A - BT, where A and B are fixed. In that case 
equation (19) is really 


T* = 


P x + Af 
Bf 


Hence T is a function only of Pp f, and the material. Then if t} 

and H, defined by equations (9) and (l6), respectively, replace y and 
T, the moment-urn equation becomes 


'Vi 


, . (^m\ „(! . Tt 




( 20 ) 


while the appropriate energy equation is , again, equation (17b ) . 

Now the forcing term (the right side of eq. (20)) should be consid- 
ered. According to its definition (eq. (16) ) H is an increasing func- 
tion of T. Hence, at least for constant viscosity, the forcing term is 
an increasing function of H. If p and k are constant, then the 
forcing function is linear in H and, for this case, it has been shown 
(refs. 1 and 6) that two solutions occur under certain conditions. In 
the present situation things are not that simple, and the results depend 
on how p and k vary. For many liquids the conductivity varies only 
moderately over a fairly wide range of temperature, while the viscosity 
may change severalfold. Two cases, water and liquid sodium, are illus- 
trated in tables I and II. In the present discussion the variation of 
conductivity is neglected. The viscosity can be written to good approx- 
imation as 


H = 


s 


T + T c 


( 21 ) 


where 
and 


T a are constants (for water s = 0.36 (centipoise) (°C) , 


s and 

Cb 

= 20° C if T is in °C). This expression is compared in tables 

I and II with experimental values. It is to be expected that 

(T + T a ) >0 in the range where the fluid remains a liquid. Under these 

circumstances, equation (16) yields 


H = T(T + 2T a ) 


( 22 ) 
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and equation (20) becomes 


U T}71 



(23) 


From equation (17b), it is seen that, if u^ is large, H must 
vary more or less parabolieally upward across the channel. Then it fol- 
lows (eq. (23)) that, when H - H(T*) has a large magnitude, the forcing 
term increases only slowly with H, while for small H - H(T*), the 
forcing term is linear in H. The latter reduces the problem to the 
usual free- convection situation (ref. l), while the former (large 
H - H(T )) approaches the usual Poiseuille case, wherein the forcing 
term is constant. This circumstance at least restricts the range of 
flow parameters for which two solutions, as found in reference 1, can 
exist. 


Stated more explicitly, if (T - T Wo )/(T Vq + T a ) remains small, the 

viscosity is essentially constant and the system becomes that solved in 
reference 1. 


Had the variation of conductivity been allowed for, a small modifi- 
cation of the foregoing argument might occur. If the conductivity drops 
an the temperature rises, the forcing term would move toward a more lin- 
ear variation with H. 

In general, the solution of the system given by equations (23), 

(l7b) , and (8) is not simple. However, after two limiting cases are dis- 
cussed, the. general case can be described. The first such case is that 
of small frictional heating; the second is for small temperature varia- 
tions, and therefore viscosity can be considered constant. 


Small Frictional Heating 

First consider the case of small frictional heating. The formula- 
tion involving q and H (eqs. (23) and (l7b) ) is not convenient for 
this case. Hence, consider equations (18), (7d) , and (8). Hie viscosity 
is defined by equation (2l), and the conductivity is assumed constant. 
When the frictional heating is small so that the right side of equation 
(7d) is negligible, one obtains 

T ■ ^0 + 1 ( V - V 


(24) 
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Then equation (18) can he integrated to give 


U » — 


p[3fd 2 

24s 


- \ 0 ) z [3(f) 4 + Cta 2 + 


Cl2“lf 2 - “3>(t) - 2 c ^3@1 


(25) 


where 


% - * 




31 ’ V - % 


tXo = 


T w 0 + T a 

^1 - \ Q 


12a^cx2 + 4(02 + 2<x-j_) + 3 
a 3 = Cl + 2a ? ) 


(26) 


and T* is defined by equation (l9) . From these, the wall shear is 

(pV y=0 = P f P d ^! " %> (^A 2 ) 

(nuy) Y= & = pfpdC^ - T^q) (03/12 - c^L - 1/2) 


(27) 


The net mass flow is 

.d 


/ 

•'n 




1 + 


5(1 + 2o 1 )(l + 6a 2 + 6cc|) 


1 + 


2<x 2 


(28) 


These results can be compared with those for constant fluid properties. 
(This case is given on p. 10 of ref. 1.) In such a case, the temperature 
is again given by equation (24) . The velocity distribution, wall shear, 
and mass flow are given, respectively, by 


“ - - V> Cl t -3) [I (M(I- ♦ -x)] 


(29) 
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(M’Uy ) Y=o - P f 1 3ta -(' I w 1 - J(l + 3 ccq.)/6J 
^Vx=d = pfpdCT^ - T W() ) [-(2 + 30^)76] 


f 


pu ffl = 


p 2 d 5 fg 

48s 


( T W! " T w 0 


)*[(l + 2cx 1 )(l + 



(30) 

(31) 


In these last equations the viscosity has "been taken as that correspond- 
ing to the average temperature] that is. 


_ 2s 

^ -^Wj_ + ^Wq + ^a 


(32) 


A comparison of the results obtained for the cases of constant and 
variable properties is given in figure 1. The fluid is liquid sodium, 
and the temperatures of the two walls (100° and 900° C) differ enough 
that the viscosity varies by a factor of about 4. In spite of this there 
is no significant difference between the results for the two cases. That 
is, the effect of variable viscosity is unimportant even though the tem- 
perature variations are large. 


It Is interesting to observe -the case when T* = (T w ^ + T VJ .^)/2, the 
average fluid temperature. Then op = -l/2, 0,3 = -1, and the shear Is 


not only the same at each wall, but is the same in the constant- and 
variable- viscosity cases. However, the velocity profiles differ slight- 
ly, and, while the mass flow is zero for the constant-property case, it 
Is not for variable properties. 


Small Temperature Variation 

In the case where the frictional heating is considered but where the 
temperature variations are small, the fluid properties can be considered 
constant. This case has been solved by Ostrach in some detail (refs. 1 
and 6), by machine methods. However, another method of getting the same 
results is now presented that has the advantage of giving the parametric 
dependence simply. The same method Is applied later to the general case 
(large temperature variations). However, the justification of the pro- 
cedure is most convincingly displayed by comparison with the foremen- 
tioned machine solutions. 

In this case of small temperature 'differences, it Is again conven- 
ient to work from equations (18), (7d), and (8). If the temperature 


46bl 
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changes are small enough, that (T - T W q)/(T W q + T & ) is small everywhere, 
the viscosity and conductivity can he considered constant. 

The following dimens ionless quantities are defined: 


where 




t = (T - T*) — 


P 2 p 2 f 2 d 4 _ K(T - T*) 


16(5^ - arj 


2 t - 1 


(ppfd 2 ) (T^ - T*) 


is the parameter defined in reference 1. In terms of the new variables, 
equations (18) , (7d) , and (8) become 

^T]T) = -T (35) 

>1 - - a n ( 56 > 

U(±l) = 0 

t(-1) m K/l6 \ 

t(+1) = h (^r = ^ 16 ’ say I (37) 

An approximate solution of this system can be found by iteration. The 
velocity profiles are usually parabolic. Thus, suppose 

= A(1 - T} 2 ) (38) 

where A is an undetermined parameter. Then equation (36) yields, sub- 
ject to the boundary conditions, 

t = -fi? [(m + 1) + i)(m - 1)J + A " ■ ^ 1 (39) 
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If this is put into equation (35) , the resulting velocity distribution 
is 



This process could be continued, assuming convergence, but is stopped 
at this point. If equations (39) and (40) are put into equation (36) 
and the result is integrated across the channel, a quadratic equation 
for the unknown parameter A 2 follows. This is 


17,325 / K(m + 1)\ 


• - 
1 m 

S08 y- 210 j 

+ 103,424 

4(m + l) 2 


(41) 

Real solutions exist if 


210 


a/_Z07_ Q _ 
\/ 2,640 


< K(m + l) < 


210 


m 


(m + l ) 2 


-i , . / 707 r m 

2, 640 [ 4 " (m+ x) 2 


(42) 


The boundaries defined by equation (42) are plotted in figure 2 for 
K > 0. Some li mi ting values found in reference 6 by machine methods are 
shown for comparison. Agreement is excellent. 


The two solutions can be examined generally in the following manner. 
If K(m + l) is moderate (say, in the range 0 to 40), then equation (41) 
yields, approximately, the two results 


or 


A 2 « jo. 016 K(m + l) 


42 


(41a) 


The mass flow and heat transfer to the walls are, respectively, 

,d 

n,i flY « ^5_ 

24 


/ pu ay « K(SL + .1) 

djjj 24 


or 


Pfd 

« ==- 16 
pfd 


(43) 
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A( ^ k m T Y^ = \i ^Y^Y=0 “ ^%^Y=d] 

- T*) 

~ 3 ^ 0.016 K(m + 1 ) 


(44) 


or 


256V(Tw n - **) 


3Kd 


40 


The first (smaller) solution is one in which there is negligible friction- 
al heating and, hence, heat transfer, while the second is quite the oppo- 
site. The second case is, as is pointed out in reference 4, one of regen- 
erative heating. There is a large amount of heat transfer to the walls 
(eq. (44) ) j this heat is supplied by frictional heating of the fluid 
occasioned by large mass flow (eq. (48)) and the resultant high shear. 
Notice that even when T^ — = T , so that the problem is an homoge- 

neous one (K = 0 in eq. (37), although K/ (T^ - T*) in eq. (33) is not 

zero) , this second solution does not vanish. For that matter, the second 
solution is virtually independent of K and m, provided K and m take 
on moderate values. 


The range of validity of these results is limited by the condition 
that (T - ®wq)/( t wo + he small so that the viscosity variation is 

negligible. In the case of the second solution (large A^) the maximum 
of T occurs near the center of the channel, and thus equations (38), 
(44) , and (46a) yield 


•'max - -"wq 

V> + T a 


£00 % - ** 
K ^wq + ^a 


If, for example, K = 10, then T v ^ - T is limited to a few degrees. 

This implies also that the wall temperatures must be virtually equal. On 
the other hand, for applications of the first solution (small A^) , the 
only restriction is that (T w ^ - T W ^)/(T W ^ + T a ) be small. However, if 

the viscosity is that for some suitably defined average, the error due to 
larger temperature variations should be unimportant. This conjecture is 
based on the example discussed after equation (32) . 


Before giving a numerical example, it Is worthwhile to examine the 
order of magnitude of the numbers one obtains in physical problems for 
the present case of Rma.ll temperature differences. From equations (33), 
(34), (40), and (41a), one has, very roughly. 


ft 
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u(t] = 0) 


m + 1 
16 


L 

'*(*» 0 - T *>' 

V- 

1 


or 


25 


/k(T W0 - T*) 


|iK 


(Small solution) 


(Large solution) 




T(ti = 0) - T V() ~ (^f 1 ) (^ Q - T*) (Small solution) 
or 

- lj (T Wq - T*) (Large solution) 

For example, for water at 0° C, if K = 25 , _ T V q - T* = 0.25° C, and 
m = 1, the numbers for the various cases are 


(45) 



Small 

Large 

u(tj = 0), ft/sec 
T(q = 0) - 'E Wq , °C 

25 

0 

200 

2 


This case corresponds to a channel width of 14.6 inches. A more general 
idea of the orders of the numbers involved in the second solution can be 
obtained as follows. From equation (33) 


T 


u 2 k 

_ T # M- T 


(46) 


However, from equations (39), (40), and (41a) the maximums of U and t 
for the second solution usually occur near q = 0 and have the respective 


values of about 6^ and 14. Hence, for the second solution, 


uf_\ 

- 


3k/n 


(46a) 


Actually there is no combination of K and m such that equations (41), 
(39), and (40) yield U^ BX / < 1* For water, 3k/ p, is of the order 

10^(ft/sec) 2 /°C, and hence if T - T* is 1° C, the velocity maximum is 
100 feet per second. For liquid sodium, 3k/ |i is about 10^(ft/sec)^/°C; 
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thus a temperature difference of only 1° C corresponds to a maximum ve- 
locity of 3,000 feet per second. It therefore appears that, if the ve- 
locity is to be kept moderate to maintain laminar flow, the temperature 
variation must be small. Hence, the assumption of constant fluid prop- 
erties is a good one. 

To compare the present method of calculating the velocity and tem- 
perature profiles with the more exact solutions obtained by machine meth- 
ods (ref. l), a single example is shown in figure 3. Remember that this 
is a constant- fluid-property situation. The example is that of water 
flowing in a channel 14.7 inches wide and for which T w ^ - T* = l/l0° C. 

The wall temperatures are 20.0° and 20.1° C. This leads to K = 10. The 
agreement is within 10 percent for the second solution and, not surpris- 
ingle, virtually exact for the smaller one. 

Before ending the discussion of exactly fully developed flows, it 
should be observed that the iterative procedure used to get solutions 
here can be applied, in the other cases considered in references 1, 4, and 
6, namely those involving wall temperature gradients and heat sources in 
the fluid. 

As is stated earlier, the present iterative procedure can be applied 
directly to the original problem wherein the frictional heating is con- 
sidered and large temperature changes are contemplated. The solution for 
such a flow is given in appendix C for the case of equal wall temperatures. 
The only difference from the case just discussed is that some of the in- 
tegrals are rather involved and the equation for the amplitude is more 
complicated. However, as is observed previously, if the velocities are 
to be kept moderate in the second solution, the temperature variation will 
be negligible. 


Approximate Solutions 

The results thus far presented have the beauty of being exact within 
the limitations of fully developed flow. However, several cases arise in 
which such a limited view is unacceptable. The simplest such case is the 
flow of a gas with a pressure gradient, the ordinary Poiseuille flow. An- 
other case of some interest is that involving a wall temperature gradient. 
The extent to which such flows can be considered fully developed is ex- 
amined in what follows. 

We seek solutions of equations (l) to (6). Again, a long channel is 
assumed and end effects are neglected. In such a case, any gradients in 
the X-direction (in the flow direction) must be small. Hence, write the 
variables as follows: 
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P = P |l + eg(x) + . . .] 


p = p 

P 0 (y) + 8pi(x,y) + . . . 

T = T 

T 0 (y) + ST 1 (x,y) + . . . 

- ■m 

p = p 

M-qCy) + 8tii(x,y) + • • • 

“ — 

k = k 

ko(y) + &k 1 (x,y) + • . . 

u = u 

u Q (y) + 8u 1 (x,y) + . . ._ 

a 

v-E 

: u jo + Sv-^x^) + . . .J 


x = X/L 
y = Y/d 


(47) 


where d/L is small, d being the wall spacing and L being a length 
of flow, as yet undefined. The other two parameters, e and 8 are 
small but unrelated at this time. For a gas 8 must be at least as 
large as e for the state equation (eq. (5)) to make sense, while for 
a liquid the number 8 is determined by the temperature boundary condi- 
tions. The barred quantities are, except for u, given parameters chosen 
so that p Q , kQ, Pq, and so forth, are of unit order. The value of u 

is initially unknown because there is no characteristic velocity for in- 
ternal flows of this kind. 


Had the term g(x) in the pressure (the first of eqs. (47)) been 
considered as a function of y also, then added terms would be intro- 
duced because of the y-momentum equation (eq. (3)). However, these are, 
analogously to the usual boundary-layer analysis, of higher order than 
■what is retained in the other equations of motion. 

If equations (47) are put Into equations (l) , (2), and (4) and only 
the dominant terms of each kind are retained, there follows 

p_(u_ + v n ) + v n + p, u_ = 0 

i,x l,y' K 0,y 1 K l,x 0 

HiT [ p 0^ u l,x + Vl^y)] + | F -- £ (X) + pf ( p 0 + & Pl)| = p 



( 49 ) 
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SpC^uT 


fpoCuQ^ljX + v l T 0,y^| + L |^ u l, x * Y ljy|] d 2 V 


S(koTi, y + ilT 0 ,y)yl + R^“0,y)l 

-I d L -I 


Finally the state equations (eqs. (5)) yield 

P = pRT 


P 0 T 0 1 


PqT-l + p-jTq = (e/&)g(x) 


P 0 - 1 = -PT(T 0 - 1) 
p x = -pot-l 


1 


(Gas) 


(Liquid) 


(50) 


(51) 


(52) 


Now observe that the terms in the braces in equation (49) describe 
the driving forces of the flow and must therefore be of the same order 
of magnitude as the viscous terms (which contain the highest deriva- 
tives) . The following four cases will be considered in turn: 


Case 

Body force 

Large temperature variation 
(Tq variable) 

I 

Yes 

YeB 

n 

No 

Yes 

in 

Yes 

No 

IV 

No 

No 


Case I - Tq Variable, Free Convection 


For case I, terms of order 5 can be neglected as compared to cor- 
responding zero-order terms. Thus, p^_ can be set equal to zero in equa- 
tion (49), as can k-j_ and T-j_ in equation (50) . Assuming p 0 and 
g x (x) to be of uni t order, the^ body force and pressure gradients must be 
of similar size so that 



( 53 ) 
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Actually, it is only necessary that — £.< pf in order to have the body 

force matter. However, it can be assumed that equation (53) holds, and 
e can be discarded later if it is small enough. 


Then, using equations (51) and (52) , the driving term in equation 
(49) is 


g x U) + pf(po) “Pf 1 + S v (x) - 


L 6 x 


1 


T 0 - 


Tr 


l" 


(Gas) 


or 


, p#[ [iliW - (* 0 - i)‘ 

. pT 




(Liquid) 


(54) 


For a mixed flow, g x (x) = 0(l) for a gas or g x (x) = -1 + 0(pT) for a 

liquid. A pure free- convection flow might arbitrarily be defined as one 
for which g x (x) = -1, but, for convenience, any flow involving body 
forces is henceforth referred to as a free-convection flow. 


Definitions of the barred reference values in equations (47) are all 
straightforward, except that for_ u, and a selection of values can read- 
ily be made a priori. However, u must be chosen such that uq is of 
unit order, and there is no way of knowing ahead of time how big the 
flow will be. Hence, for the moment, let us beg the question and define 
simply 

C 2 p - PfpT (55) 


where Cg is unknown and, for a gas, (3T = 1. It is shown later that 
C 2 is a number in the range 10 to 50. Then equations (48), (49), and 
(50) can be written 



Po( u l,x -* 

• v l,y> + p 0,y v l + Pl^ " 0 ( 56 ) 

j 

i T n - 1 

r , 1 ) 8u 2 c 2 r 1 

^0 u 0,y)y + C 2 j 

T° - 

1 + g x (x)||. |p 0 (» 0 Ui ;X +v lU0)y )J 

\ 


(Gas) 

^O^py + °2 | 

o* 

I 

H 

i 

[ + s x «] / ST | = ^ [p 0 (”o”i, x + n"o,p] 


(Liquid) J 
(57) 
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tic v 5 u 2 C 2 


fLpT 


CyPT 

pe v Su^Co ( T 1 

(k 0 T 0,y)y - — — =- | L P ° C “° Tl »* + T l T 0,y>J 


| PQv'^l.x + v l T 0,yjj + 
fp 0 (uoTi jX + Vi"o, y ) - “0 | 8*00 | 


mT ^O^y) 


(Gas) 


f t*,* + v l T 0,y)_ } - ^ (■V’O.y) (11,luia) 


(58) 


The system whose solution is sought is given by equations (56) , 

(57), (58), and (51) or (52) plus the boundary conditions (eq. (6)). 
Assuming that g x (x) is given, the unknown dependent variables remain- 
ing are seven in number, Uq, u-j_, Vj_, Pq, p-j_, Tq, and Tj_. With five 

equations and seven dependent variables, some further restriction must 
be made. The difficulty arises mainly in the terms u-j_ and v-j_. For 
a gas it appears p^_, uj_, and v-j_ are of the same order as Tj_ and 

hence that the inertia term in equation (57) is of the same size as the 
convection term in equation (58). Hence, it follows that, in order for 
the solution to be determined, both of these terms must be negligible. 

The same result follows for liquids, although perhaps not so obviously. 
In this case (eq. (5l)), p^ is very small, of the order of PT, where 

P is the volumetric expansion coefficient. However, for liquids the 
viscosity is a very strong function of temperature, and thus = 0(Ti). 

Hence, there is no particular reason to assume that u-j_ and v^_ are 
not the same size as T]_. Accordingly, if these terms matter, the pres- 
ent formulation is useless. 


For these terms to be negligible two courses are open. One is to 
have everything x- independent as in the exact solutions described ear- 
lier. The second is that the parameter CjjSu^/fLpT be small. This is 
not simple. For example, if equation (55) is used and a channel 1 inch 
wide and having a characteristic length L of 10 feet is assumed, then 
under standard conditions and gravitational acceleration 

C 2 5u 2 /fLpf = 66,000 5/C 2 
=6,400 S/0 2 


(For water ) ) 
(For air) ) 


( 59 ) 
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Consider Cg = 1. For the air case, B must be as large as e, which is 
about l/3,000 here. For water B can be chosen by an applied wall tem- 
perature gradient. In either case it seems difficult to make CgB^/fLpT 
small, in fact, unless Cg is a large number. 

If this question is ignored for the moment, it can be observed for 
a gas that, if 8 = e (and 8 must be as large as e), the coefficients 
of the inertia or convection terms and of the dissipation terms are 

Cg8u 2 /fLpT = C 2 eu 2 /fL = Cgpu 2 /P 

and 

12 - p- 

kT c p 

(where Pr is a Prandtl number) which are both essentially the squares 
of a Mach number and are the same size if Cg = 1. Hence, for the gas 

case, at least, the frictional heating must also be negligible because, 
as is shown later, Cg > 1. 

For a liquid it can be seen, by trying some cases, that the fric- 
tional heating must again be very small provided S/e is not virtually 
zero . 

Now to return to the question of a value for Cg. Suppose that equa- 
tion (55) is used together with its subsequent consequences. Then, for 
the sake of an example, let g x (x) = -1 and assume a gas. The differen- 
tial equations become 

C^O u O,yV = c 2^ 1 " T o )/ T 0 
fro T 0,y)y = 0 

The boundary conditions are 

UqCo) = UpCl) = 0 

T 0 C°) - T »o 



ToCU = 


TCQ* 
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To make the point about orders of magnitude, consider that pq and 
kg are proportional to the temperature (i.e. , M-q = 3£q = Tq) . Then the 

equations can be solved very readily. Maximum speeds have been computed 
for several cases and are shown in the following table: 


Tw o 

^Wq_ 

UO 

C 2 

^ max 

1/4 

1/2 

1 

1 

3/2 

2 

0.12 

.01 

.03 


Them j uq! max is of urLit order lf c 2 is > sa &> 40 * Klus for the condi- 
tions cited in connection with equation (59) , 

C 2 5u 2 /fLpT = 1,6505 (For water) 


or 


= 1606 (For air) 


which, particularly for the air case, can readily be made small. In a 
case (such as eqs. (35) to (37)) where two solutions occur, the same con- 
clusion about C 2 follows, although the argument is rather tortuous . 


Finally, then, the inertia terms must be negligible for the flow to 
be f ully developed. In that case the thermal convection and frictional 
dissipation are negligible. The equations for fully developed free- 
convection flow in a channel are (eq. (55) is assumed to apply with 
C 2 = 1) 


r \ 1 - T o 

vl\)' u 0,y/y ~ 

^W^y = ° 


(Gas) 


(61a) 


1 + g x 

fcW’o.yb - 1 - T o + -^=- 

(*0*0, Py = 0 


(Liquid) 


(61b) 
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Cg ss 1 is used here and also in equation (55) because C2 was intro- 
duced only as an aid in determining what matters in the equations of 
motion. With these small terms eliminated, Cg can be dropped. 

It should be remembered that these equations correspond to cases 
where the temperature variations across the channel can be large. The 

solutions are valid provided only that Su^g^/fLpT is small. For con- 
sistency, of course, g x (x) must be constant. 


The system is readily solved by first integrating the energy equa- 
tion. The result in parametric form is 



where R is the right side of the first of equations (61a) or (61b). 

This then is the solution of a free -convection flow wherein large 
temperature variations across the channel are admitted. The effect of 
longitudinal wall temperature variations would presumably be allowed for 
by considering that these were local profiles, by a sort of strip theory. 
Results of the kind in reference 4 would apply to case III, discussed 
later on, wherein small temperature variations are assumed throughout. 


Case II - Tq Variable, Forced Convection 

This is a flow in which the body force is considered to be negligi- 
ble. In such a case equation (53) no longer applies, but one assumes 
that f = 0 and that e is given. Without loss, take g x (x) = 1. 
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Equation (55) is replaced by 


Pe _ pu 
L " d 2 

Then equations (6l) and (62) with their accompanying conditions hold, but 
with the right sides of the first of equations (61a) and (61b) replaced 
■by unity. The solution is that given in equations (62) but with R = 1. 

The analysis thus far given applies to the case where large temper- 
ature differences are allowed (Tq ^ l) . If only small differences are 
permitted, a somewhat different formulation results. 



Case III - Tq = 1, Free Convection 

If Tq = 1, then pQ = Pq = 1. Also, equation (53) still applies, 
and hence the driving term in equation (49) is (compare with eq. (54)) 


Pe 

L g x 


(x) + pf (1 + 5 Pl ) = pfS [-T, + -|S + 1 

_ r i + g-x-tx) 

= pf8pT -T-, + z^r— 

|_ x 5f3T 

Then, for a free-convection flow, it is required that 

g x (x) = -1 + O(SPT) 


(Gas) 
(Liquid) 


l (64) 


(65) 


Under such a condition the driving and viscous forces are of the same 
order provided (compare with eq. (55)) 


pf5pT = pu/d 2 


( 66 ) 


If these last equations and equations (5l) and (52) are used, equations 
(48) , (49) , and (50) become 


u l,x + v l,y " -Pl^UQ 


(67) 


+ T , - aM - 1 + ^ (x) , 

66T fLpT 


“0,jy + T 1 - sr~ - ^ ("ov + v i"o,y ) 


( 68 ) 
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^ - Si t 1 “HJ ' « ^ C69) 

where the terms involving e/6 do not appear if the fluid is a liquid. 


In distinction to the result for the case of large temperature vari- 
ations , in equation (68) , 

“2 o 

— = 66,0005^ (For water) 

fLjBT 

or 

= 6,4006 2 (For air) 

for the circumstances of equation (59). In the case for air, if 6 = e, 
this is a very small number and should, indeed, "be negligible. The con- 
vection term in equation (69) is of order 6 larger, as is the friction- 
al heating term, and should be retained. 

Then equations (68) and (69) are 


”0,KT + 0 ' - T 1 + 1 + - 

+ S x (x) 
83T 

(70) 

M 

Tl ' yy ’ S + [rtf ( x + f^) 

T l,x ' SxJ "oj 

(71) 

- 


where again the terms involving e/6 in equation (70) and g x in equa- 
tion (71) appear only for a gas. 


Now to determine what form the temperature variations can take. It 

was assumed at the outset that Uq is independent of x. For this case 

it is readily shown that the most general forms allowable for the temper- 
ature and pressure gradient are 

Tq = a-jx + a 2 + T 2 (y) (72) 

g x (x) = -1 + &PT £(e/8 + aq)x + a g j (73) 

Then equations (68) and (69) become 


%yy ~ t 2 


(70a) 
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I 


jxu 2 


2,yy " ^ i-^,y + u o 


1 + a, — rr 
^ fLp 


+ *§J 




(71a) 


where the term. (iu 2 /&kT { Uq [l] } disappears for a liquid. These equations 

require small temperature variations but admit substantial mass-flow 
rates . 


Then, in terms of Uq and Tg(y) , equations (70a) and (71a) are a 
pair of ordinary differential equations and are nonlinear only if the 
frictional heating is imp ortant. Solutions in the linear case are quite 
simple (ref. 4) and in the nonlinear case can be found by the iterative 
method given earlier following equation (37). Ostrach discusses this 
system extensively in reference 4, where among other things, some machine 
solutions are given. 


Case IV - T 0 = 1, Forced Convection 

For case XV equation (63) still applies, and equations (70) and (71) 
are replaced by 


"o^yy “ 


g x (x) 


(74) 


■ W K 


SpTc 


+ 


Fe 


1 (l + =M + 


} 


(75) 


The temperature variation must have the form 

T i = ( a i + a 2^ x + *2(y) 


while g x (x) is a constant. These equations are easily solved because 
the dependent variables are separated. The result is, of course, the 
familiar Poiseuille flow. 


Other Gas Flows 

The flows discussed thus far have included only those cases wherein 
the velocity is essentially independent of distance along the channel. 

If this restriction is lifted, the problem becomes vastly more complicated. 
Therefore, only one class of solutions is examined here. These solutions 
are ones in which the independent variables are separable. Such a flow 
is out of the question for liquids unless it can be assumed that the tem- 
perature is a function of y only. This is because of the form of the 
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state equation (eq. (5b)). It can be shown that this limitation to T 
independent of x leaves only the fully developed cases discussed 
previously. 


First recall that, to have a fully developed flow in any sense, the 
channel must be very long and the dependence on x must be much weaker 
than that on y. Then equations (l) to (5) can be appro ximat ed as 


(pu) x + (pv)- 

Y ~ 0 

(76) 

p(uu^ + VUy) + P X = 

-Pf + (HUy)^. 

(77) 

II 

o 


(78) 

pc v (uT x 4- vTy) + PCux + ' 

Vy) = (kTyJy + 

(79) 

P = pRT 

(Gas) 

(80a) 

p = p |l - p(T - T)] 

(Liquid) 

(80b) 

These equations can be derived formally in the same manner as they are 
derived for external boundary layers. The only difference is that the 
Beynolds number of boundary-layer analysis is replaced by a ratio (l/d) 2 , 
where L is a characteristic flow length and d is the channel width. 

If the fluid is a gas (eq. (80a)) and the body force is negligible, 
and the viscosity and thermal conductivity each vary with the temperature 

? 

ii 

ZL 

9 


(81) 

IV 

II 

IT* 



then the permissible separated forms are 


CD 

W" 
' — * 

ii 

v > 

5* 

'B' 

6X/L 

or e ' 


v/v-l(Y) = (X/L) 9-1 

or e fflc/li 


p/ Pl (Y) = (x/L) ®( 2 ^" 1 ) +1 

or e (2t-l)SX/L 

K (82) 

t/t-^y) = (x/l) 26 

or e 20X/L 


P/P-L = (X/L) e ( 2 £ +1 ) +1 

or e (2J+l)0X/L 


where 9 is an arbitrary constant. 

✓ 
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It is interesting, though perhaps irrelevant, that exactly the same 
variations of free-stream velocity are allowed for similar solutions of 
the external-hound ary-layer equations (ref. 7). 

The exponential form in equations (82) is valid for the complete 
Navier-Stokes equations, while the other form depends critically on the 
assumption of a very long channel. 


Several other somewhat unrelated comments about this result are per- 
haps in order. First, for a liquid the requirement previously stated, 
that Sr/SfX = 0, leads to the condition that nothing varies with X. 

This case has alr eady been examined. Second, for a gas , if the viscosity 
and conductivity do not have the same variations with temperature (eqs. 
(8l)), o nly the trivial X- independent separation results. The 
X- independent solution corresponds in equations (82) to the exponential 
variation with 6=0 and was discussed starting with equations (7). 
Finally, if the body force is important (-pf in eq. (77)), the forms 
given in equations (82) apply, but only with 6=0 (exponential and un- 
interesting) or 6 = l/2 (power of X). 

The forms given in equations (82) have two other properties of in- 
terest. The through- flow Mach number, which is proportional to u/ -/T, 
is Independent of X. Also, -unless there is flow through the channel 
walls, all the solutions except 6=0 (exponential) and 6 = -l/2£ 

(power of X) must be flows with zero net mass flow. This is because 
the mass flow is 


or 


f pu dr 
u 


(X/L) 2 ^ 1 f p 1 (Y)u 1 (Y)dY 

e 2£6X/L /’ d p (y) u ( y)<3Y 

J 0 


( 83 ) 


which must not vary with X unless there is flow through the walls. 
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If equations (82) are put into equations (76) to (80a) the result 
is 

^±A) PiUi + ( Pl v l)y = o 

p l (l U 1 + v 1 u 1,y) + + P 1 “ W U 1,Y>Y - Pl f 


P 1,Y - 0 


>-(84) 


C v p l (ir U 1 T 1 + V 1 T 1 ,y') + P 1 (l U 1 + v 1,y) " (^1®1 ,y)y + Miu| #y 

P 1 = P 1 RT 1 

where A = 0 if the exponential variation is used and A = 1 if the 
variation is as a power of x. Also* the term p^f can appear only if, 

as already mentioned, 6=0 (exponential) or 0 = l/2 (power of x) . k-, 

and have the obvious definitions. x 

Now a new space variable, q, is introduced from equation (lO). Then 
equations (84) become 


BP 1 (2^0 + A) 
2RLT^(-1) 


U 1 + = 0 


(85) 


R 


{[^^] *}* - - ¥ 
1 - 2 s) E - Ej p l}“l T l + T* 

4bT^(-l) ]__f 

m -P / 


( 86 ) 


B 


— T i C*L,im + f 


p i = pi rt i 


(87) 

( 88 ) 
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Observe that equations (85) , (86), and (87) form, a set of three ordinary 
differential equations for three variables (uq, Tq, and pqvq) . One 
curious feature of the system is that it is of fifth order and there are, 
in general, six boundary conditions to be applied (one each to uq, vq, 
and Tq at each vail). Hence some restriction must be placed on the 
combinations of boundary conditions for given values of the parameters, 
particularly Pq and 0. Actually this is no different from vhat hap- 
pens in analysis of the external laminar boundary layer. In that case, 
hovever, the bo undar y condition on v (or v-^) at the outside of the 
boundary layer is not satisfied, nor is there any particular reason for 
it to be. For channel flovs such an omission is probably not allovable. 
In general, this means simply that there must be a flov through at least 
one w all , and this flow cannot be prescribed if similarity is to be 
maintained. 

If, on the other hand, one examines the case where 2£0 + A = 0 
(eq. (85)), it is permissible to set v-^ = 0. Then the system is of 

fourth order with four boundary conditions. The condition 2£0 + A = 0 
corresponds exactly to the two cases mentioned in connection with equa- 
tion (83) . These are the only ones in which v-^ can vanish at the 

walls. Of the two cases defined by 2£0 + X = 0, namely, X = 0, 6=0 
(exponential x- variation) and X = 1, £6 = -l/2 (x®- variation) , the first 
has been solved earlier. In the second case, for linear viscosity- 
temperature variation, equations (86) and (87) become (-20 = £ = X = l) 


u^ + RTq = 


- - gs 


Tin 


Jet l 1 - n?) 


^(KTP - + f 


where the boundary conditions are 

Uq(H) = 0 

Tq(-l) = Constant 
T^(+l) = Another constant 


(89) 


(90) 


(91) 


The solution of this system is more difficult to obtain than is the solu- 
tion of equations (35) to (37), although a similar procedure can be 
followed. 
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If the flow is fairly slow, with the Mach number limited to, say, 
a few tenths, the profiles are, to good accuracy 


T 1 - 





U 1 aLCT^+l) + Txi-1)) 


s 


Pi d 



- I X (-1) \ nfl . 

+ ^T^+l) + T^-l) 1 12 




where use has been made of equation (10) to determine B. It may be ob- 
served that the distance L can be defined by 


L = - 



which follows directly from the last of equations (82). If the wall tem- 
peratures are equal, the velocity profile is exactly the familiar Pois- 
ieulle one. 


A final remark: It can be seen that the system described by equa- 

tions (90) and (91) will probably admit pairs of solutions just as the 
free-convection flow of equations (35) to (37) does. 


CONCLUDING REMARKS 

When fully developed channel flows are considered, the cases that 
can be solved exactly are very limited. For a gas a constant pressure 
is required, or at least one which does not vary in a streamwise direc- 
tion. This case is analogous to Couette flow in that no approximations 
need be made in arriving at a relatively simple mathematical problem. 

In the case of a liquid, one can solve the exact case of constant pres- 
sure gradient in the streamwise direction. For both the gas and the 
liquid, the wall temperatures must be constant. In the gas case nothing 
astonishing happens. However, in the liquid case a surprising result 
arises. There appear to be (except for certain singular cases) either 
two or no solutions for the flow. This result, which has been discussed 
extensively by Ostrach, has one solution for which frictional heating is 
negligibly small. The second is one in which the frictional heating is 
large, and thus the temperature is raised and the buoyancy effect is in- 
creased. In the present report an approximate analytic solution of this 
problem is given. The results agree very well with Ostrach’ s machine 
calculations.. Although an analysis is given for the case of variable 
viscosity, It turns out that for the cases of interest, wherein the fluid 
velocity is kept within reason, the temperature variations are small and 
there is no reason to consider variable viscosity or conductivity. 
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These so-called exact solutions, particularly in the case of a gas, 
do not cover all the flows of interest. Hence, consideration is given 
to cases in which there are streamwise temperature and pressure gradients 
but in which the flow velocity is virtually independent of distance along 
the channel. For a gas the mere presence of a pressure gradient requires 
a temperature gradient, while for a liquid the presence or absence of a 
streamwise temperature variation is governed by the wall temperature 
conditions . 

In these cases one of two situations occurs. If the temperature 
variation across the channel is of the order of the temperature level, 
then in order that "channel flow" be maintained, the convection terms in 
the energy equation must be negligible. This implies that the mean flow 
Mach number is small and also that the frictional heating is negligible. 
For such circumstances the equations are separated and can readily be in- 
tegrated for any case of interest. Only one solution exists. 

On the other hand, if the temperature variations are small, more 
complicated effects occur. This situation of very small temperature 
changes admits very large flow velocities (see the discussion following 
eqs. (46)). Then both the frictional heating and the thermal convection 
effects can be significant. In such cases (when the frictional heating 
matters) two solutions can occur. These flows qualify as quasi- 
incompressible in that the only place where compressibility effects mat- 
ter is in the buoyancy term in the streamwise momentum equation. 

The forced flow perhaps deserves an added remark. When the temper- 
ature varies only slightly across the channel, the velocity profile must 
be the usual parabolic one. When the variation in temperature is large, 
the profiles can still be found in closed form but are more complicated. 

If streamwise variations of velocity are allowed, the flow is more 
complex. A description is given of the circumstances under which the in- 
dependent variables are separable. These forms can yield new results 
only for gas flows and show that the streamwise variation must be either 
as a power of x (streamwise coordinate) or exponential ly with x. With 
two exceptions, only one of which admits x- variations , these flows re- 
quire that the body force be negligible. The exponential causes apply to 
the full Navier-Stokes equations, while the other ones require an expan- 
sion of the equations of motion in terms of the width- to-length ratio of 
the channel. For all these cases, the streamwise Mach number is inde- 
pendent of x. All but one of these possible flows lead to difficulties 
with boundary conditions and require a flow through the walls. The lone 
exception has streamwise velocity proportional to l/ -y/x. For small flow 
Mach numbers the solution is similar to that for Poisieulle flow but al- 
lows for temperature variations across the channel. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 

Cleveland, Ohio, June 5, 1958 
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APPENDIX A 
SYMBOLS 

parameters defined In eqs.. (41) and (C5) 
parameters defined In eqs. (9) 
parameter defined in eqs. (10) 
parameter defined in eq. (C9) 
specific heats 
vail spacing 

body forces in X- and Y-directions, respectively, considered 
positive in the minus X- and Y-directions 

pressure perturbation, eq. (47) 

temperature functions, eqs. (l6) or (22) and (Cl) 

parameter defined in eq. (34) 

thermal conductivity 

characteristic length, eq. (47) 

parameter defined In eq. (37) 

pressure 

gas constant, eq. (5a) 
parameters defined in eq. (21) 
temperature 

reference temperature defined in eq. (19) 

dimensionless velocities in eqs. (33) and (Cl) , respectively 

velocity components in X- and Y-directions 

Cartesian coordinates, X being in the main flow direction 
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I 

8 j€ 

n 


p 

p 


x/L and Y/d, respectively 

parameters defined in eqs. (26) 

volumetric expansion coefficient (see eq. (5b)) 

small parameters introduced in eqs. (47) 

dimensionless distance across channel in eq_. (10) 

parameters defined in eqs. (81) and (62) and after eq. (84) 

viscosity- 

density 


t dimensionless temperature difference in eq. (33) 

Subscripts : 

m mean value corresponding to average of wall temperatures 

Vq wall conditions at Y = 0 

Wj_ w a l l conditions at Y * 4 

X,Y,x,y,Ti partial derivative with respect to that variable 
0,1 zero-order and first-order solution in eq. (47) 

Superscripts : 

(o) , (l) first two approximations 

reference values in eqs. (5b) or (47) 


bars 
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APPENDIX B 


A GENERALIZATION OF EQUATIONS (7) 

If a transverse tody force f-[_ and a distribution of heat sources 
pQ, where Q is a constant, is included, equations (7) become 


v = 0 

(pu y )y = Pf 

Py = -P^ 


■> 


o 

(kI Y ) Y = -pu£ - pQ 


(Bl) 


These equations can be solved in exactly the same manne r as equations (7). 
The results are," for linear variation of viscosity and conductivity with 
temperature, f ^ 0, and Q = 0, 



TQQ* 
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where 


a = BCi/283^ is defined hy 


a 


3g _ _s5_ /£!oR\ 2 r, ff 
d f l bdf i \^2af| J If 


e" CT ) 2 (l/3 - </* 2 ) + ^ ( e 2a 



and 


, for f-j_ = 0 and Q ^ 0, 


P = Constant 


u.^UZ.D 




r (B3) 


% + ^1 f Tw l" Tw 0 ) a /aB“fP 

2 2 j - b ( ^ 

^0 




-hc- » - '■> 

8hRT£ 


where 


1 Wa 


240 Id 3* 


f 2-,3\ 
Pm fd ] fd 

^m k m / ■ C m 


r B^m \ / Pm^ d 


V ia A/ A^o, 


■ 1 


If Q, is positive, there is only one real root of this equation, that 
root being such that 


If Q is sufficiently negative, there can be three real positive roots. 

Other solutions can readily be obtained for the case where neither 
f^_ nor Q, vanishes or where other distributions of heat sources occur. 
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APPENDIX C 


EXACT SOLUTION FOR LIQUID 

The problem at hand is to solve equations (27) and (17b), subject 
to equations (8) . Define 


Hi = 


krf + 2sH(T) /t + T 


a 


3sT| + 2sH(T*) V t * + T £ 


Ui = u 


2s 


u-i/2s/k 


y 


kTf + 2sH* T* + T 

a< a. 


(ci) 


In terms of these variables the problem is 


= _ / B^fsV2i7k \/ _ _1_ 
^ ( 4(T* + T a ) 11 


H 1>iTO = " 

^(± 1 ) = 0 

Hi(±l) = [(^ Q + T a )/(T* + T a )] 


(C2) 

(C3) 

(C4) 


where, to keep the problem from getting out of hand, equal wall temper- 
atures are assumed. For convenience, define 

A 1 = [(% + Ta>/( T * + *a>] 

(C5) 


B^pgfsV2s/k 
4(T* + T a ) 


It should be remembered that the parameter B, and hence Ag, is as yet 
undetermined (recall eq. (lO)). 


The solution is found in exactly the same manner as the solution of 
equations (35), (36), and (37). Thus assume 


u(°) = A(l - T) 2 ) 
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Then equations (C3) and (C4) yield 


H-L = A x + A 2 (l - n 4 )/3 


If this is put into equation (C2) , that expression becomes 


u- 


(1) _ 
•i>Trn 



(C6) 


(C7). 


Tf this is integrated and the appropriate boundary conditions are satis- 
fied, the result is 




(i) 


= + V^/A 



dTi 


'o 


_ / ,1 __d ! L_ + f 1 JU2L 

J 0 A 


where 


C 4 = 1 + 3A.J./A 2 


(C8) 


(C9) 


If equations (C8) and (C6) are put back into equation (C3) and the result 
is integrated across the channel, the result is 


8A 2 = 2 _ 2-V5 
3A| 3 A 


f 1 (1- n 2 )dn 1 + j. r( r 

Jq V c4 - a2 Jo 


_dTL 


Vc 4 - 


dri (CIO) 


n 4 


Before using this equation to determine A, something has to be done about 
B, which is as yet undetermined. If equations (2l) , (Cl), and (C6) are 
used, the second of equations (10) yields 


t3A(T* + T a ) 

V^sB 



(Cll) 
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If this is put into the first of equations (C5) , eliminating B, and 
the result is put into equation (CIO), the defining equation for A 
is, finally. 



I 

This can he solved to very good accuracy by 



(C12) 


(C13) 


The second root of equation (C12) is never a, significant one. This equa- 
tion can readily be solved for A in terms of C (which is itself a 
function of A; see eq. (C9)) and the parameter 



K 

V^(V% - 1) 


where K is defined in equation (34) . As in the case of constant prop- 
erties, two possible values for A are again found. 

The various integrals appearing in equations (C8) and (C13) can be 
evaluated as follows (ref. 8) : 
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where F(cp,l/-/2) and E(cp ,l/-/2) are the respective elliptic integrals 
of the first and second kind of amplitude cp and modulus 1 / -^2- 
of these integrals can he approximated as follows 


= -4? (tj + t] 5 /10C 4 ) 


I C 4 - Tl 4 C 


which is correct to 1 percent if C > 1.3, and 



which is correct to 1 percent if C > 1.0. 
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TABLE I. - VISCOSITY AND CONDUCTIVITY 
OF LIQUID SODIUM 


Temper- 

ature, 

T, 

°C 

Conductivity, 

K 

watts/(cm) (°C) 

Viscosity, 

U, 

centipoises 

Experi- 
ment 
(ref. 9} 

1117 
1170 + T 

Experi- 
ment 
(ref. 9) 

Eq. (21) 
(a) 

200 

0.815 

0.815 

0.450 

0.450 

300 

.759 

.760 

.345 

.350 

400 

.712 

.712 

.284 

.286 

500 

.668 

.668 

.243 

.242 

600 

.627 

.636 

.210 

.210 

700 

.590 

.598 

.186 

.185 

800 

.547 

.567 

.165 

.166 

900 



.150 

.150 


= 150° C; s = 158(centipoises) (°C) . 


TABLE II. - VISCOSITY OF WATER 
UNDER ATMOSPHERIC CONDITIONS 


Temper- 

ature, 

T, 

°C 

Viscosity, 

V-, 

centipoises 

Experi- 
ment 
(ref. 9) 

Eq. (21) 
(a) 

0 

1.79 

1.79 

20 

1.01 

.90 

40 

. 66 

.60 

60 

.48 

.45 

80 

.36 

.36 

100 

.28 

.30 


®T a = 20° Cj 

s = 36 (centipoises ) (°C) . 






Figure 1. - Effect of property variation on velocity profiles. 
100° C; , 900° C; T*, 0° C; liquid sodium. 
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Figure 2. - Regions of existence of solutions of equations (35) to (37). 
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Figure 5. - Velocity and temperature profiles, m, 2j k, 10j 
T^ - T*, l/lO° C j water. 
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